Energy Level Diagrams for Black Hole Orbits 



o 
o 

(N 

> 

O 

(N 



i , 



(N 
> 
in 

On 

in 
o 

On 
O 



X 



Janna Levin 

Department of Physics and Astronomy, Barnard College of Columbia University, 3009 Broadway, New York, NY 10027 and 
Institute for Strings, Cosmology and Astroparticle Physics (ISCAP), Columbia University, New York, NY 10027 

janna@astro.columbia.edu 

A spinning black hole with a much smaller black hole companion forms a fundamental gravitational 
system, like a colossal classical analog to an atom. In an appealing if imperfect analogy to atomic 
physics, this gravitational atom can be understood through a discrete spectrum of periodic orbits. 
Exploiting a correspondence between the set of periodic orbits and the set of rational numbers, we 
are able to construct periodic tables of orbits and energy level diagrams of the accessible states 
around black holes. We also present a closed form expression for the rational q, thereby quantifying 
zoom-whirl behavior in terms of spin, energy, and angular momentum. The black hole atom is not 
just a theoretical construct, but corresponds to extant astrophysical systems detectable by future 
gravitational wave observatories. 



Bare black holes are as perfect as fundamental parti- 
cles. As Chandrasekhar said, "The black holes of nature 
are the most perfect macroscopic objects there are in the 
universe [H". A black hole with a given mass and spin 
is indistinguishable from every other black hole with the 
same mass and spin. Likewise, a supcrmassivc black hole 
with a much smaller black hole companion forms a kind 
of macroscopic, classical atom, reminiscent of the hydro- 
gren atom. In analogy with atomic physics, the orbits 
around a given black hole can be completely described 
by a periodic table Q - a table of periodic orbits ordered 
in ascending energy from the stable circular orbits (the 
ground-like state) up to the last bound orbits (the energy 
of ionization). Further, the energy levels of the periodic 
orbits around a black hole are, formally speaking, dis- 
crete. In the spirit of this atomic analogy [? ] , we exploit 
energy level diagrams as another valuable representation 
of the dynamics around a black hole nucleus. 

The energy level diagrams are based on an infinite or- 
dered set of periodic orbits around a Kerr black hole [1] . 
Fig. Q] shows a periodic table of retrograde orbits. At 
first glance, the set of idealized, closed orbits might seem 
improbably special. Yet the periodic set is fundamental 
to an understanding of the entire dynamics. As we now 
elaborate, all orbits are near a periodic orbit, and there- 
fore this seemingly special set is able to encode the entire 
orbital dynamics. 



PERIODIC ORBITS AND THE RATIONAL 
NUMBERS 



Our atomic analogy will stem from a mapping between 
black hole orbits and rational numbers. For the sake of 
argument, wc consider orbits confined to the equatorial 
plane. Any equatorial orbit around a black hole has two 



fundamental frequencies 1 - the libration in the radial co- 
ordinate, ui r , and the rotation in the angular coordinate, 
u> v . To be clear, by ui r we mean the frequency of ra- 
dial oscillations, cu r = 2n/T r , where T r is the time it 
takes to move from one apastron to another. By lu v , 
we mean the average angular frequency over a radial os- 
cillation, uj v = J 1 " 1 Jq 1 ' (<kp/dt)dt. A generic orbit will 
have frequencies that are irrationally related and will 
consequently never close. Instead, it precesses around 
the central black hole nucleus, never returning exactly to 
its starting point. In contrast, a special set of orbits - 
the periodic orbits - will close. For this to be the case, 
the two fundamental frequencies must be commcsuratc, 
that is, rationally related. We define the rational number 
q associated with each periodic orbit through 
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Since the set of irrationals is a larger infinity than the set 
of rationals, the rationals are a set of measure zero in the 
space of all possible orbits. 

However, the set of periodic orbits are dense on the 
number line. Any irrational can be arbitrarily well ap- 
proximated by a nearby rational. A good way to find the 
rational approximant is to use continued fractional ex- 
pansions. Let a be some irrational number and let 1 + a 
be the irrational ratio of frequencies for some generic or- 
bit. Such a trajectory precesses around the central black 
hole without ever closing. The irrational can be repre- 
sented through continued fractions as 
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ai 



(2) 



0.2 + - 



where w is the integer part of a, ai is the integer part 
of the reciprocal of a — w, is the integer part of the 



1 Non-equatorial orbits have three fundamental frequencies and 
therefore are designated by two rationals. 



2 




FIG. 1: A periodic table for retrograde equatorial orbits 
around a black hole with spin parameter a = —0.9. The 
angular momentum is L = 4.56056. All entries with z < 4 and 
< w < 2 are included. Energy increases from top to bottom 
and then left to right, as does the rational q — w + v/z. The 
columns are w-bands 0, 1, 2 from left to right. The minimum 
is q ^ 0, as discussed in £11111 and tjlVI Although we have 
drawn 3 columns, there are in principle an infinite number as 
w — > oo. 

reciprocal of the remainder of a\ etc. through an infinite 
list of integers a n . A given rational approximant is ob- 
tained by truncating the continued fraction at some finite 
n. For instance, if we truncate the continued fractional 
expansion at n = 3, then 



z aia 2 a 3 + oi + a 3 

In this way, any irrational can be arbitrarily well approx- 
imated by rationals; significantly then, any aperiodic or- 
bit (with an irrational frequency ratio) can be arbitrarily 




FIG. 2: Left: A generic retrograde orbit around a central 
black hole with spin parameter a — —0.995. The angular 
momentum L is halfway between the ibco (innermost bound 
circular orbit) and the isco (innermost stable circular orbit). 
This orbit has an energy very near the 3-leaf orbit on the 
right. 



well approximated by some periodic orbit (with a ratio- 
nal frequency ratio). 

To illustrate, the generic aperiodic orbit on the left of 
Fig. [2 if allowed to run indefinitely, would fill an an- 
nulus in the plane, never closing. This is a generic ex- 
ample of zoom-whirl behavior. (For other discussions of 
zoom-whirl behavior see 0.) Yet this same orbit has a 
a w q = 2 + 42/125, and so is very close to a periodic 
orbit that whirls 2 times, has 125 leaves but skips to the 
42nd next leaf in the pattern each successive time it hits 
apastron. Also note that the orbit is still close to, but not 
as well approximated by, the much lower-leaf orbit char- 
acterized by the rational q = 2 + 1/3, drawn on the right 
of Fig.O Indeed, the generic orbit can be thought of as a 
precession of this fundamental three-leaf pattern by the 
amount 27r(42/125 — 1/3) = 27r/375 - so less than one 
degree each radial cycle. In fact, in general, all high- 
leaf orbits can be understood as precessions of lower-leaf 
orbits 0. 

For emphasis, every orbit that is measured observa- 
tionally or generated numerically is necessarily truncated 
by a rational approximant and is therefore, strictly speak- 
ing, periodic, albeit possibly with an extremely long pe- 
riod. Consequently the orbital dynamics around the cen- 
tral black hole is entirely defined by the periodic skele- 
ton. Furthermore, high-leaf periodics can be understood 
as precessions of lower-leaf periodics. 

Although we restrict to two-dimensional equatorial 
motion here, generic three-dimensional non-equatorial 
orbits can also be identified with rational numbers. In 
principle, two rationals are needed to characterize a com- 
pletely closed non-equatorial orbit, the q we are using as 
well as a qg = ujg/oj r . However, it should be mentioned 
that we do not need to require fully closed motion, as 
shown in Refs. [? ? ] in the context of comparable 
mass binaries. The useful insight is the simple observa- 
tion that every orbit lies in its own orbital plane, even 
if it does not lie in the equatorial plane. Periodic tables 
of orbital plane motion were shown to taxonomize orbits 
with one rational, <?$ where $ is the angle swept out in 
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the orbital plane, just as the periodic tables of equatorial 
orbits are taxonomized by the q above. The entire orbital 
plane then prccesscs to trace out fully three-dimensional 
motion, although not necessarily fully periodic motion. 
For the rest of this paper, we will restrict to equatorial 
motion for transparency. 



II. q AND ITS TOPOLOGICAL SIGNIFICANCE 

We can relate the triplet of integers (z, w, v) in the 
definition of the rational 



■ v/z 



(4) 



to geometric and topological features of the orbit. Con- 
sider Fig. [1] again. The periodic orbits are assembled in 
a table of increasing energy - as is the chemical periodic 
table. In our tables, energy increases from top to bottom 
and left to right. Note that the q of the periodic table 
increases monotonically with energy. 

Each q also has a geometric interpretation. Consider 
its definition 



v u, 
q = w H — = 

z U> T 



•P _ 1 = ^f. _ 1 
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(5) 



where Aip = J Tr (dip/dt)dt is the equatorial angle accu- 
mulated in one radial cycle from apastron to apastron. 
By this definition, we see that q is the amount an orbit 
precesses beyond the closed ellilpse. Therefore Kepler's 
elliptical orbits have q = since there is no precession. 
Keplerian orbits are all periodic. They are also relativis- 
tically excluded. Relativistic orbits around black holes 
all overshoot their starting point, which is to say they all 
precess, a discovery that originates with Mercury's famed 
precession, and all relativistic orbits have a q > 0. 

Note that all of the orbits in the second column of the 
periodic table of Fig. [T] whirl an additional 2tt in a nearly 
circular loop around the nucleus before moving out to 
the next apastron. The entries in the third column whirl 
twice. The integer part of q from equation ([5|), denoted 
by w, is the number of whirls. The periodic tables are 
ordered in terms of «;-bands, which form energy bands. 
In fact, unlike the chemical periodic table, the black hole 
periodic table in Fig. [1] has infinitely many columns, one 
for each possible integer number of whirls. 

Within each energy band, each w-band, the entries are 
z-leaf orbits that whirl w times between each apastron. 
Closed orbits can skip leaves each radial period. For 
example, the third and fifth entries of column 2 in the 
periodic tables are both 3-leaf clovers but they are not 
identical. The third entry in column 2 moves to the next 
apastron in the 3-leaf pattern as indicated in bold. How- 
ever, the fifth entry in column 2 moves to the second 
apastron in the 3-leaf pattern as indicated in bold. Ac- 
cording to equation ([5]), the q of the former closed orbit 
is q = w + v/z = 1 + 1/3 while the q of the latter is 
q = w + v/z = l + 2/3. So v labels the order in which 
the leaves arc hit. 




84 
0.915 



= 2 + 2/3 
= 2 + 3/5 



= 2 + 2/5 
= 2+1/3 



FIG. 3: Energy Level diagrams for equatorial orbits with L = 
2 around a central black hole with spin parameter a — 0.995. 
All z < 5 for all 1 < to < 12 are included although the q = 1 1 
and q = lj entries are missing from the w = 1 band since 
they are not accessible. The figure on the bottom shows a 
detail of the top figure, indicating the rationals within the 
w = 2 band. 



The rational number identifying a given periodic or- 
bit therefore encodes significant topological and geomet- 
ric information. It identifies the number of leaves, the 
number of whirls, and the order in which the leaves are 
traced. 



III. ENERGY LEVEL DIAGRAMS 

Since the rationals are discrete, it follows that the ener- 
gies of the periodic orbits arc likewise discrete. Further- 
more, as discovered in Ref. [2J, q is monotonic with energy 
(for a given angular momentum L around a given black 
hole). This latter observation allows us to represent the 
information in a periodic table in a related energy level 
diagram as shown in Fig. [3] 

Energy level diagrams provide a valuable complement 
to the periodic tables, effectively and concisely summa- 
rizing information. For instance, Fig. [3] incorporates all 
z < 5 orbits for 1 < w < 12, while the corresponding pe- 
riodic table of Fig. [T] only includes z < 4 for < w < 2. 

The energy level diagram illustrates the stacking of 
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energy levels for all orbits into bulk uj-bands. The frac- 
tional part of q, represented by v/z, has an intriguing 
distribution within a given band that roughly repeats in 
a scaled manner throughout all w-bands. 

If we were to draw every one of the infinite g's allowed, 
they would, despite being discrete, fill a solid wall of color 
for each w-band. Still, even though there are always an 
infinite number of entries in any periodic table, high-z 
orbits can be approximated as precessions arounds low- 
z clovers. In practice, we can restrict to a finite set of 
discrete rationals. 

Furthermore, not all rationals are accessible for a given 
angular momentum, L, around a given black hole. The 
lower limit is set by the q of the ground state - equiva- 
lently, the q of the stable circular orbit - and the upper 
limit is set by the q of the highest-energy bound orbit. 
Both upper and lower limits depend on L and the spin 
of the central black hole. 

Although initially counter-intuitive, the lower limit is 
not simply 0. As in Ref. Q, we take the q for the ground 
state, that is, for a given stable circular orbit, to be 
q c — > Lu v /uj r — 1 in the limit of eccentricity going to zero. 
This minimum q is not zero. For the energy level diagram 
of Fig. [3l the bound on q is roughly 1 + 1/3 < q < oo. 
Importantly, this means that some low q orbits are not ac- 
cessible in the strong-field regime. Namely, Mercury-type 
orbits, which as precessions of Keplerian ellipses would 
correspond to rationals < q <C 1, are forbidden in the 
strong-field regime and all allowed states will look like 
precessions of multi-leaf orbits. We will also write down 
an explicit analytic expression for q c in ijlVI 

As the angular momentum increases both the lower 
and upper bounds will eventually decrease toward zero, 
so that, as we know from solar system dynamics, all orbits 
in the Newtonian approximation are tight precessions of 
the ellipse. 



where an overdot denotes differentiation with respect to 
the particle's proper time r and 



R(r) = P 2 - A {r 2 + (L - aE) 2 } 
P(r) = E(r 2 +a 



aL 



A 



r 2 - 2r + a 2 



(8) 



where the mass of the black hole is set to M = 1. 
We then have 



1 



{aE-L) + fP 



n 



dr - 1 



(9) 



where we take the plus sign in f indicating we are inte- 
grating from periastron, r p , out to apastron, r a . 

To solve this, we write both R and A in terms of their 
roots. R is a fourth-order equation in r. In the equatorial 
plane, one of the roots is always and R can be written 

R{r) = (E 2 - l)r(r - r )(r - r p )(r - r a ) (10) 

while we can write A in terms of the inner and outer 
horizons r± = 1 ± y/l — a' 2 : 

A= (r-r + )(r-r-) . (11) 

We want to solve 



I r r " dr 1 

q = ( a E -L) --= + -a 2 (aE - L) 

TT J rp VR X 

1 r» r 2 dr 

+ -aE / = 

7T J rv AVR 



which wc rewrite as 



dr 

aVr 

•1 , 

(12) 



IV. CLOSED FORM EXPRESSION FOR q(a,E,L) 

We can always numerically integrate the Kerr equa- 
tions to find q. However, it is useful to have a closed form 
expression purely in terms of (a^E.L). Such a closed 
form should be advantageous when considering the evo- 
lution of q under the dissipative effects of gravitational 
radiation [? ] . Starting with the definition 



q = — f —dr — 1 
2tt.L r 



(G) 



we can perform this integral to generate an expression in 
terms of elliptical integrals. We will take the further step 
of finding the limits of integration in terms of (a, E, L) 
so that our final expression will automatically deliver q 
given (a, E, L) as the only input. 

We rewrite @ the geodesic equations for Kerr equato- 
rial orbits [1] as 



r 2 r = ±VR 



r 2 ip 



}_dR 
~K~dL 



(7) 
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Defining 



--(aE - L)X X + -a 2 (aE - L)l 2 + -aEX 3 

TT TT TT 

dr 
1 



dr 

aVr 



AVi? 



a(r) 
(3(r) 
7(r) 



the solutions are 



r{r p - rp) 
r P (r - r a ) 
r(r p - r a ) 
r P (r - r a ) 
(r - r a )(r a - r p ) 
(r a - ro)(r - r v ) 



(13) 



(14) 
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Mr) 

Mr a) 



y/(E2 - l)r p (r - r a ) 
2 

- l)r p {r^~) 



EllipticF 
EllipticK [P(r )\ 





1 




ArcSin 









(15) 



with Ii(r p ) = 0; The second integral is 



Ti = 



-2 



W(E 2 -l)r P (r a -r ) 



( EllipticF 



ArcSin 
ror+ 



(r - r_)(r_ - r+) 
Tor- 



a(r)- 1 / 2 ] ,a(r ) 
EllipticPi 
EllipticPi 



(r - r+)(r+ -r_)" 

This gives non-zero contributions at both limits of integration; Finally, we have 

2Q> - r ) 



'airy 1 ' 2 


,a(r a ) 


'air)' 1 ' 2 


,a(r a ) 



Is =" 



y/(E* - l)r p (r - r~) 



(r - r_)(r p - r )(r - r+) 



EllipticF 



ArcS 



(ro - r_)(r_ - r p )(r_ - r+) 



EllipticPi 



(»'o - r+)(r+ - r p )(r+ - r_) 



EllipticPi 



7 (r_), ArcSin 
7 (r_|_), ArcSin 



,j9(ro) 
1 



1 



,0(n>) 
,/?(r ) 



(16) 



At r = r a , j(r a ) = 1 and an EllipticF becomes an EllipticK, while at r = r p , 7 (r p ) = oo and EllipticF and EllipticPi 
become zero. So evaluating this at the upper and lower limits gives 



Mr a) 



2(r p - r ) 



- l)r p (r - r a ) 



(ro - r_)(r p - r )(r - r+) 



(ro - r_)(r_ - r p )(r_ - r+) 



EllipticK [/3(r )] 



EllipticPi [ 7 (r_),/?(r )] 



We use these to evaluate 
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(r -r + )(r+ -r p )(r + -r_) 
1 



EllipticPi [ 7 (r + ),/3(r )] 



1 



9 = (oE - L)Ii(r ) + -« 2 (a£ - L) [I 2 (r ) - X 2 (r p )] + ~a£X 3 (r a ) - 1 

7T 7T 7T 



(17) 
(18) 



We are not quite done since this expresses q(a, E, L, ro, r p , r a ). To have q as a function of (a, L) only, we have to 
find ro, r p , r a as functions of (a, E, L). To do so we expand the polynomial R(r) = [E 2 — l)r(r — ro)(r— r p )(r — r a ) and 
equate to the definition of R(r) in Eq. ([8]), matching up coefficients in powers of r and finding a system of equations 
for r , r p , r a . Since r = is always a root, this is equivalent to a 3rd order equation in r and cubic equations have a 
generic solution. We want the roots of 



r) = 







D = 







A = 


(E 2 — 


1) 


B = 


2 




C = 


Aa 2 - 


L 2 


D = 


2(L- 


aEf 



(19) 
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The nonzero roots in ascending order are 



where 



r __ _ B_ + (1 -^ + 3,C) (F + ^_^ {F+ ^ 



F = -2B 3 + 9 ABC - 21 A 2 D 

G = F 2 - 4(B 2 - 3AC) 3 . (21) 

Using Eqs. (|2"D|) in Eq. (fT5)) , gives the rational g as a function of (a, 25, L), which is what we were after. 

To be clear, the function q(a, E, L) of Eq. (|18[) can in principle be any positive real number. The periodic orbits 
correspond to rational q. All other orbits can be arbitrarily well approximated by a rational q as discussed previously. 
As a check, inputting the (a, E, L) of the retrograde orbit on the right of Fig. [2] returns a q = 2.3333... = 2|, as it 
should. 



While Eq. (fT5|) gives the g for a generic equatorial orbit, 
it is interesting to write out the expression for the lower 
limit of stable circular orbits explicitly. Circular orbits 
have a radius r c = r a = r p and this leads to a great sim- 
plification in the expression since then (3(ro) = 7(r+) = 
7(r_) = and EllipticK[0] = EllipticPi[0, 0] = tt/2. Us- 
ing these relations gives 



1 



v/(£ 2 -l)r c (r -r c ) 



-(aE-L) 



aE 



A(r c ) 



The root r$ corresponding to a given r c can be found 
from the cubic expressions above evaluated at the energy 
and angular momenta of circular orbits, originally found 
in Ref. [? 1 



E = ±- 



3/2 „ 1/2 , 

r c — 2r c ± a 

3/4, / 3/2 o 1/2 , o 



L = ±- 



r 2 T 2ar, 



V2 , „2 



3/4 / 3/2 „ 1/2 , „ 

?V \/r c -3r c ± 2a 



(22) 
(23) 



where the ± refers to progradc and retrograde orbits re- 
spectively. 

Alternatively, and perhaps more easily, we can write 
the roots as functions of (a, E, L) by expanding, 

R(r)/r = 

= [E 2 - l)(r - r )(r - r c f 

= (E 2 - 1) (r 3 - r 2 (r + 2r c ) + r{r\ + 2r r c ) - r r 2 ) 

to find the r ^ roots when r a = r p are degenerate at 
the value r c . Equating coefficients with those of Eq. (fTI 



10 r 




2 4 6 8 10 12 14 ' C 

FIG. 4: g c for circular orbits for spin values from left to right 
a = 0.998, 0.8, 0.5, 0, -0.5, -0.8, -0.998. Negative spin values 
are equivalent to retrograde orbits. The value of q c diverges 
at the innermost stable circular orbit (ISCO) for each value 
of the spin as expected 0. 



gives 



A = (E 2 - 1) 
-A(r + 2r c ) = 2 
A{r 2 c + 2r r c ) = Aa 2 - L 2 
-Ar a r 2 c = 2(L - aE) 2 



(24) 



and 



ro 



3(£ 2 - 1) 
2 



2r r 



x _^J) {[E2 _ l)a2 _ L2) 



(25) 



When solving the quadratic, the sign consistent with 
( r c — r o) > is chosen. Note that for all bound orbits 
E 2 - 1< 0. 
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The value of q c for different spins is shown in Fig. [4] 
In the strong-field, the ground state has a q c well above 
zero, ruling out Mercury-type precessions in favor of the 
more extreme zoom-whirl precessions. 



V. GRAVITATIONAL RADIATION 

In addition to being fundamental, the black hole atoms 
have the added benefit of being real. Black hole pairs 
formed by tidal capture in dense regions, such as globular 
clusters [? ] or the galactic nucleus [? ? ], will likely 
become bound on eccentric orbits. Light black holes that 
fall onto supermassive black holes in the galactic nucleus 
are also free to be eccentric. Supermassive black holes 
with stellar mass black hole companions are candidates 
for direct detection through gravitational radiation by 
LISA (Laser Interferometer Space Antenna). As entries 
in the periodic tables, we can identify the q of these black 
hole inspirals as a function of (a, E, L) given the results 
of this paper. This allows us to predict the relativistic 
zoom-whirl behavior from a simple formula. 

In our periodic tables and energy level diagrams we 
have neglected the dissipativc effects of gravitational ra- 
diation. As energy and angular momentum are radiated 
away in the form of gravitational waves, there will be 
transitions from the orbits of one energy level diagram 
at a given L to those of another energy level diagram at 
lower L. 

The rate of change of q can be expressed in terms of 
the change of energy and angular momentum as 



dq 
~dt 



dq dE 
dE~dt 



dq dL 
dL~dt 



(26) 



This can be evaluated by explicitly taking the derivatives 
of (fT8|) . being sure to treat rn, r p , r a as functions of (E, L) 
through Eqs. (|2"D|) . Accurate calculations of dE/dt and 
dL/dt for zoom- whirl orbits require detailed numerical 
computations @, @, 0]- I n conjuction with such compu- 
tations, our closed form expression for q would enable 
us to find possible resonances during inspiral Q, which 
occur during episodes of dq/ dt w 0. As seen in Fig. 
dq/dE > and dq/dL < 0. Consequently it is possible 
that the condition Aq m will be satisfiable: 



dq 
~dE' 



AE 



dq 
dL' 



AL 



(27) 



during inspiral and an orbit will noticeably hang near a 
given q for a time. While a detailed study will require a 
full numerical treatment of adiabatic inspiral, we mention 



that in the simulations of Ref . Q , it can be seen by eye 
that the orbit will occasionally hang on a very nearly 
closed orbit as it sweeps through to merger. 

In full numerical relativity we have seen q change 
rapidly for comparable mass black hole pairs [? ]. By 
contrast, an extreme- mass-ratio system can spend hun- 
dreds to thousands of orbits in the strong-field regime 
where periodic tables are most illuminating. The small 

30 „ — „ — , — „ — ,— — „ — , — „ — ,— — r - 



25 - 
20 - 
<S> 15 - 
10 - 

5 - 



0.82 0.83 0.84 

Energy 




1.85 1.90 1.95 2.00 2.05 2.10 

Angular Momentum 

FIG. 5: a = 0.995 Top: q versus E for L half way between 
the ibco and the isco, L — 1.81058. Bottom: q versus L for 
E = 0.852079. 



black hole will cascade through a sequence of orbits, each 
of which looks like a precession of a low-leaf clover. As 
it inspirals into the black hole nucleus, it should leave an 
imprint of its transit through the periodic table in the 
gravitational radiation. 
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